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Abstract 


A weakly nonlinear theory is developed to study the interaction of TS waves and Dean 
vortices in curved channel flow. The predictions obtained from the theory agree well with 
results obtained from direct numerical simulations of curved channel flow, especially for 
low amplitude disturbances. At low Reynolds numbers the wave interaction is generally 
stabilizing to both disturbances, though as the Reynolds number increases, many linearly 
unstable TS waves are further destabilized by the presence of Dean vortices. 
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I. Introduction 


The study of incompressible curved channel flow can illuminate some important issues deal- 
ing with laminar- turbulent transition. The streamwise curvature of the flow induces an 
instability resulting in longitudinal vortices (hereafter referred to as Dean vortices [1]). For 
the case in which the gap between the walls is small relative to the radius of curvature, Reid 
[2] described two approximate methods which both showed that the vortices first become 
unstable when the Dean number (De = 2Re^/ d*/r*, where Re = U*h*/v is the Reynolds 
number, U * is the bulk velocity, d * is the channel width, h* = d * / 2 is the channel half-width, 
rj is the inner wall radius, and v is the kinematic viscosity) exceeds a value of about 36. 
These vortical structures have been found experimentally by Ligrani and Niver [3]. The 
linear and nonlinear properties of these axisymmetric vortices as well as their transitions 
into two distinct types of wavy vortices have been well documented in the direct numeri- 
cal simulations and weakly nonlinear studies of Finlay, Keller, and Ferziger [4]. For more 
mildly curved channels than they studied, the minimum Reynolds number for instability of 
the Dean vortices is comparable with that for two-dimensional (2D) Tollmien-Schlichting 
(TS) waves. In this regime, Gibson and Cook [5] showed that oblique waves are never the 
dominant linear disturbance. Hence a study of the nonlinear interactions between 2D TS 
waves and the streamwise oriented Dean vortices can be meaningfully conducted. 

Daudpota, Hall and Zang [6] (hereafter referred to as DHZ) developed a weakly nonlinear 
interaction theory to study the interaction of Dean vortices and TS waves in curved chan- 
nels. They employed a multiple scale version of the Stuart [7] and Watson [8] approach to 
derive two coupled Landau equations for the perturbation amplitudes of Dean vortices and 
TS waves. A comparison of their theory’s predictions with the results of direct numerical 
simulation [9] suggests that their results are in error with respect to the influence of the 
TS wave on the Dean vortex. The resolution of this discrepancy is the major motivation 
for undertaking this work. In section II. we present a slightly different formulation of a 
weakly nonlinear interaction theory and describe how the resulting equations are solved. 
The current approach is based on the work by Herbert [10, 11]. We compare some aspects 
of this approach with that used by DHZ. In section III. we feport some results of the inter- 
action theory while in section IV. we compare specific cases with results of direct numerical 
simulation. Finally in section V. we draw conclusions. 
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II. Mathematical Formulation 


A. The Basic Equations 


The incompressible Navier-Stokes equations in cylindrical coordinates ( r*,6,z *) are written 
as: 
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and p* and u* represent the constant density and kinematic viscosity respectively. The 
asterisks indicate dimensional quantities. The geometry of the problem of azimuthal flow 
between infinite concentric walls of outer radius r* and inner radius r* is illustrated in Figure 
1. The wall boundary conditions require that: 


u r * = Us * = u; = 0, at r* = r*, r* 0 . 
A solution to the equations gives ( U r *, U$*, U 2 * f P*) as: 
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Here U* is the bulk velocity and 
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The quantity r* c = (r* + r*)/2 is the centerline radius and the ratio r) = r*/r* describes the 
channel curvature. 


We nondimensionalize all spatial coordinates with the channel half-width, velocities with 
the bulk velocity, and pressure with p*U* 2 . The temporal scale is h*/U*. The Reynolds 

number is defined as 

Re = Vh*/v*. 


The nondimensional equations are easily obtained from the dimensional ones by replacing all 
starred quantities with their corresponding nondimensional unstarred ones and noting that 
p = 1, h = 1, and v = 1/Re. The expression 1 — 7/ appears often below and is denoted by A. 
When A = A c = 2.179 X 10 -5 the minimum Reynolds number for instability of TS waves is 
the same as that for Dean vortices (see Gibson and Cook[5]). 


We follow closely the perturbation method introduced by Herbert [10, 11], extending it 
to the case of interacting disturbances. The steady solution is perturbed such that 
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where it r , itg, u z , and p are perturbation quantities. Substituting the perturbations into the 
nondimensionalized Navier-Stokes equations, subtracting the steady flow component, and 
rearranging gives: 


where: 
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The entire flow is assumed to be periodic in the azimuthal (streamwise) and axial (spanwise) 
directions; the solution evolves in time. In addition the total pressure gradient is kept 
constant. 
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B. Fourier Expansions 


The linearized disturbance equations are obtained by setting the right hand sides of Equa- 
tions 13-16 equal to zero. A solution to the equations for a 2D TS wave can be written 

as: 
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Here a. is real and represents the azimuthal wavenumber while 7 (t) = g 0 t is a real function 
defining the phase. Any growth or decay of the TS wave is expressed in terms of an amplitude 
parameter such that: 
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and ao,o is the linear growth rate. The apparently superfluous superscripts and subscripts 
are included here for consistency with the notation used later. A convenient normalization 
for this problem is 

|«r 1,0 (re,*)| = A(t) (22) 

where r c is the channel centerline. We take 

^ 1 ’ 0> °’°(r c ) = 1, (23) 

thus fixing the amplitude and phase of the TS wave. 

In a similar manner, one can write the solution for the primary Dean vortex mode as: 
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where /? is the real axial wavenumber and (( t ) = h 0 t = 0 indicates that the linear Dean 
disturbance is steady. We also write 
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where 


( 26 ) 


<LB_ 

It - K '° B ' 

The normalization is chosen such that: 

|*«(r.,«)l = B(t), 


(27) 


or in this case, 

U, 0AI ’“(r c ) = 1. (28) 

It is conceivable that a particularly inauspicious choice of the parameters Re , A, and 
wavenumbers might lead to u r 1 - 0 - 0 ’ 0 (r c ) = 0 or u r °’°’ 1 '°(r c ) = 0 and hence the normaliza- 
tion used in either Equation 22 or 27 could fail. Such an event would lead to overflow errors 
in the numerical procedure and an alternative normalization could be used. Herbert [11] 
discusses integral normalizations which are more robust but somewhat more complicated to 
implement. In practice we have never experienced any problems with the normalizations 
in Equations 22 and 27. We attribute this to the fact that the amplitudes of the radial 
velocities of the dominant TS waves and Dean vortices do not have local minima near the 
channel center. 


When the nonlinear terms (the right hand sides of Equations 13-16) are included, the 
disturbances can interact with themselves, the mean flow, each other and all their relevant 
complex conjugates. This results in the generation of harmonics, mean flow distortions, 
and various corrections to the fundamental disturbances. Hence it is natural to expand the 
perturbation variables in a double Fourier series: 
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C. Amplitude Expansion 


The system of nonlinear partial differential equations obtained by substituting Equation 29 
into Equations 13- 16 are coupled and difficult to solve efficiently. We seek a solution by 
expanding in the amplitude parameters A(t) and B(t) about the linear solutions given by 
Equations 19 and 24. The solutions will reflect the deviations from linear behavior for finite, 
but sufficiently small amplitudes A and B. In this context, all harmonics of the primary 
disturbances are considered to be forced. Since the linear TS and Dean solutions are 0(A) 
and 0(B) respectively, their first harmonics and cross terms will be 0(A 2 ), 0(B 2 ), and 
0(AB). Higher order harmonics and cross terms will be 0(A n B l ), hence it is reasonable to 
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where all double tilde terms are 0(1) except for those with n — l — 0. Heuristically, one 
can see that the exception with the n = l — 0 term comes about from the fact that the 
lowest order mean flow distortions are generated from the product of either TS or Dean 
fundamental disturbances with their respective complex conjugates. Since the fundamental 
disturbances are 0(A) and 0(B), the product terms which generate n = / = 0 are 0(A 2 ) 
and 0(B 2 ) respectively. Substitution of the representation in Equation 30 into Equation 29 
gives: 
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Since we seek solutions which are small deviations from the linear behavior, we now expand 
the double tilde representations in sums which are the products of ascending powers of the 
amplitude functions with coefficients which are strictly functions of the radial coordinate. 
In the limit as A — > 0 and B — » 0, the solutions tend towards the linear results. Because of 
the invariance of the original equations and boundary conditions with respect to arbitrary 
translation in the streamwise and spanwise directions and the assumption of periodicity of 
the solutions in these directions, only even powers of A and B are needed [12]. Hence, 


( r ,*)\ 


/ ~ n f l 
I U 

Ug n ' l (r,t) 
u/\r,t) 

\ P n ’\r, t) j 


= EE 

m=0 A:=0 


W\ 


/ -7 

u z n ' m 'l' k (r) 

V p n ’ m ’ l ’ k ( r ) J 


A 2 m (t)B 2 k (t). 


(32) 


By substitution into Equation 31 one obtains the full representation of the perturbation 
variables: 


f U r\ 

u e 

u z 

\ P J 


oo oo oo oo 


/ u T n ' m,l ’ k (r) ’ 

n t mj t k / 


E E E E 

n= — oo m=0l= — oo k=0 

\ p n ‘ m - l - k (r) 
exp in(a9 — 7 (t)) exp il((3z — ((t)). 


; \ I 


(33) 


In order that the wall boundary conditions be satisfied at all orders of approximation, we 
require: 
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for all combinations of n, m, and k. 

The time derivatives A/ A, BjB , 7 , and ( are also expanded in terms of A 2 and B 2 such 
that: 
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Equations 35 and 36 are the coupled Landau equations which describe the growth of the 
disturbances. The purpose of the remainder of this section is to determine the Landau 
coefficients a m *. and 


D. Solution Method 


The solution expansions in Equations 31-38 are substituted into the perturbation equations 
13-16. All terms with common exponential factors are grouped together. The simplest 
nonlinear theory is obtained by considering only those terms with Fourier exponents limited 
by: 

—2<n,l<2. 

With this restriction, it is only appropriate to consider terms in the amplitude expansions 
(Equation 32) with m — 0,1 and k = 0,1. Larger values of m or k give higher order 
contributions which should not be considered at this degree of truncation. In each group 
with common exponential factors, all terms with common powers of A and B are collected. 
This tedious process of expanding the solutions, substituting into the equations and collecting 
terms with common factors of exponentials and amplitude parameters is done automatically 
using the symbolic manipulator Macsyma. The case n = m — l = k = 0 just represents the 
basic flow which has already been subtracted out. A sequence of 12 systems of equations 
remains. Each system contains four equations and four unknown expansion coefficients 
which are functions of r. A description of each of the systems is summarized in Table 1. 
The equation system number indicates the order in which we chose to solve the equations. 
To ensure that all information needed to solve any system has already been computed, all 
systems of lower order in the amplitude parameters should be solved before those of higher 
order. 

The first two equation systems are eigenproblems of the form: 
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Eq. Sys. No. 

Order 

n 

m 

1 

k 

Description 

1 

A 

i 

0 

0 

0 

Linear TS wave 

2 

B 
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0 

1 
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Linear Dean vortex 
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0 

Mean flow distortion from TS 
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B 2 

0 

0 
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1 

Mean flow distortion from Dean 
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A 2 

2 

0 
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0 

TS harmonic 

6 

B 2 

0 

0 

2 

0 

Dean harmonic 

7 

AB 

1 

0 

1 
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TS - Dean cross harmonic 
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0 
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0 

TS - Dean cross harmonic 

9 
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1 
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0 

0 

Self- correction to TS 
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1 
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Self-correction to Dean 

11 
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Dean correction to TS 
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A 2 B 
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TS correction to Dean 


Table 1: Summary of equation systems 
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For the eigenproblems, the indices n and / may only have the values 0 and 1 with the 
constraint n + l = 1. The elements of the operator L ni * are given in Appendix A. The 
solution of the eigenproblem determines the eigenvalues (we consider here only the least 
stable ones), a 0) o + i<7o,o and &o,o + i^o,o- The eigenvectors qi^.o.oC^) and qo,o,i,o(f) are known 
to within a constant factor. The normalizations in Equations 23 and 28 are used for the two 
eigenvectors. The real parts of the eigenvalues are the zeroth order Landau coefficients. 


Equation systems 3 and 4 represent the mean flow distortion caused by the TS and 
Dean disturbances respectively. For this case, n = l = 0 and the operator can be greatly 
simplified. Here we consider only the case where either m or k equals 1, while the other equals 
0. Earlier we specified that the mean pressure gradient is constant, p 0,1,0, ° = p 0,0 ’ 0 ’ 1 = 0. 
Analysis indicates that u r 0,1,0,0 = u r 0,0 ’ 0,1 = u z 0,1,0,0 = u z 0,0 ’ 0,1 = 0. The remaining equation 
can be written as: 
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An interesting restriction to the theory results from Equation 42. Consider the associated 
homogeneous problem: 
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with boundary conditions: 


= 0 at r — Vi, ?v 


Apply the transformation: 

r 2 = ( 2 / (-(2ma 0|0 + 2kb 0f o)Re) 

to the homogeneous problem and divide by (2mao,o + 2fcio,o) to obtain: 

n 1 u i 1 d<t> i ^ - o 


( 43 ) 


with 

<j> = 0 at £ = &, 

where 


£? = r\! (— (2ma 0 ,o + 2fc& 0 ,o)Re) 
Co ~ Cl (— (2ma 0> o + 2fc& 0 ,o)Re) 


Equation 43 is the Bessel equation of order 1. If the points & and correspond to zeros of 
the appropriate Bessel function, then this free eigenmode can contribute an arbitrarily large 
component to the solution of the forced problem. Moreover, a numerically ill-conditioned 
system arises whenever and are both in the neighborhood of zeros of the Bessel function. 
A sure way to prevent this problem is to require both ao,o and 6 q,o to be non-negative. This 
essentially restricts the application of the theory to those regimes in Reynolds number - 
wavenumber space where both the TS and Dean disturbances are unstable. In this regime, 
solving for t^ 0,1 ’ 0 ’ 0 and u^ 0,0 ' 0,1 is straightforward. 

The equation systems 5—12 are all of the form: 

(L n i — (27710-0,0 + |n|a 0i0 — ^5o,o T 2 A: £>o,o T |J|&o,o ^A-o,o)Nl)c ln, m ,/,fc = R-*> (44) 

where R, is the right hand side vector of the s equation system. The right hand side vectors 
consist of nonlinear combinations of the solutions obtained at lower order. The elements of 
the right hand side vectors are given explicitly in the Appendix B. Some important consid- 
erations in solving the equations are given below. (Many of these issues are addressed in 
detail by Herbert [11] for the case of a single 2D TS wave in a straight channel.) 

With the restriction that a 0 ,o and 6 0 ,o be non-negative, the homogeneous problems asso- 
ciated with Equation 44 are unlikely to have free modes. The occurrence of any free modes 
for Dean harmonics would require that (2A: + I)i>o,o be an eigenvalue for the homogeneous 
problem with wavenumber 1(3. (Note that /io,o = 0.) Such an unlikely coincidence would 
lead to obviously nonsensical results of the numerical solution. TS harmonics and crossed 
TS-Dean harmonics would be even less likely to have free modes because both real and 
imaginary parts of the eigenvalue would need to match. 
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Equation systems 5-8 for the various harmonics require numerically solving a sequence 
of linear ordinary differential equation systems. In the case of Dean disturbances, it is 
possible that the harmonic of an unstable mode is also in the unstable regime. When 
this happens, it violates the assumption that all of the higher harmonics are forced by the 
primary disturbances. While numerical results can be obtained for such cases, the results 
are meaningless since the harmonic can grow on its own and is not simply forced by the 
fundamental. 

Equation systems 9-12 define the corrections to the fundamental TS and Dean distur- 
bances and allow us to calculate the Landau coefficients a^o, &o,i, ®o,i 7 an d &i ( o respectively. 
It is in the solution of these terms that the method developed by Herbert [10, 11] differs 
significantly from the standard Stuart [7] and Watson [ 8 ] approach. 

The right hand side vectors Ft, for s = 9-12 can be rewritten as: 

R, = Ri + **Mq n|(U o (45) 

where: 


K 9 

= ai i0 
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II 

— ih 0 i 
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* 12 

II 

<=r- 
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The correspondence between the index s and the indices n, m, Z, and k is given in Table 1 . 
Note that these equations have either n = 1, l = 0, or n = 0, l = 1. No other combinations 
of n and l appear. The solution vectors q n)Tn| i ( jk for these equation systems are the sums of 
particular solution vectors, hence it is useful to write: 

q + Xn,m.l,fc( r ) ( 46 ) 

where: 

Xn,m,l,ki.*'o) ( 4*0 

Substitution of Equations 45 and 46 into 44 leads to: 

(L n j — ( 2 ma 0 ,o + n ( a o,o — *9o,o) + 2kb 0> o + K^o.o ~ *^o,o)) M) /CaXn,m,/,Jfc = « a Mq n> o,j,o (48) 
(L ni j — (2ma 0i o + n(a 0i o — *9o,o ) + 2kb 0i o + f(^o,o — ^o,o)) M) Xn.m.i.k = (49) 

Recognizing that 

(L n ,f — (n(a 0 , 0 — *50,0) + K^o.o — iA’O.o)) M) q n ,o,/,o = 0 ( 50 ) 

when n=l,/ = 0orn = 0,/ = l, one finds by inspection of Equation 48: 

- (2ma 0 , 0 + 2kb 0l0 ) X° n ,m,i,k = ^, 0 , 1,0 (51) 
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or 


Xn,m,i,fc = -qn,o,l,o/ (2 ma 0 , 0 + 2kb 0t 0 ) • ( 52 ) 

It is straightforward to solve Equation 49 to obtain x\, m ,i,k- Substituting Equation 52 into 
Equation 46 gives: 

q “ «»/ ( 2ma o,o + 2 kb o,o) ^Ln.o.t, o(r). (53) 

In order to determine the Landau coefficients k, we use the normalizations in Equations 
22 and 27. Together with the expansions in Equations 30 - 33 they require that: 
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(55) 


and 


We already used Equations 23 and 28 to set 

u r 1,0,0,0 (r c ) = 1 

and 

Ur°’ 0X °(r c ) = 1 

hence for the normalization to be correct at all orders of approximation, 

u r 1,m,0,k (r c ) = 0 for all m,k > 0 

and 

_ q for all m,k > 0. 


(56) 

(57) 


The conditions from Equations 56 and 57 can now be used to determine k 3 . Let the 
components of vector Xn,m,l,k written, as: 


Xn ,m, i,fc - 


(™r\ 

Wq 


TD, 


(58) 


V < ) 


The conditions in Equations 56 and 57 are applied in Equation 53 to require: 

w T (r c ) — k , j (2mao,o + 2A:5o,o) u T ' ’ ' ( t c ) = 0. (59) 

Since u r n '°’ l ‘°(r c ) = 1, 

k, = zv T (r c ) (2ma 0 ,o + 2fc5 0 ,o) • ( 60 ) 

Hence the Landau coefficients ac, are determined. The correction to the fundamentals can 
be found by substitution back into Equation 46. 
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E. Numerical Considerations 


All of the equation systems are solved using a Chebyshev collocation procedure [13]. The 
left-hand side operators for all of the systems (except the mean flow corrections) are similar 
and are generated with a single subroutine. The operators for the mean flow corrections need 
to ensure that only streamwise velocity perturbations are nonzero. This was most easily done 
using a separate subroutine. The right-hand side vectors were sufficiently different from each 
other that it was more convenient to use Macsyma to write individual subroutines for each of 
the equation systems. The two eigenproblems are solved first with a global solver to obtain 
the eigenvalue spectrum. The eigenvectors associated with the most unstable eigenvalue are 
obtained with a local iterative procedure. The mean flow needs to be computed using 128 
bit accuracy for the very weak curvatures that we used. We found that 65 collocation points 
across the channel (including the endpoints) produced Landau coefficients which agreed to 
at least 3 digits with those obtained with 97 collocation points. All results presented here 
were obtained using 65 points. 


F. Comparison of Approaches 

Here we compare some aspects of the current approach with that of DHZ. Since they do not 
use the normalizations in Equations 22 and 27, they do not have the conditions of Equations 
56 and 57 and hence need some alternative constraint to obtain Landau coefficients. They 
use orthogonality conditions to obtain the Landau coefficients when both Dean and TS 
disturbances are on their respective neutral stability curves. Near the neutral stability 
curves, the linear growth rates vary linearly with Reynolds number deviation off of the 
curves. Deviations from the critical curvature parameter, A c , are also assumed to affect the 
linear growth rates linearly. They find that the effect on the coefficients of the nonlinear 
terms in the Landau equations are independent of deviations away from the neutral Reynolds 
numbers and critical curvature to the order considered. 

While the two approaches involve the same type of expansion, the different ways of solving 
for the Landau coefficients makes each appealing for obtaining certain types of information. 
Below we discuss some of the restrictions and advantages of each method. 

Both theories assume that the fundamental disturbances are the major structures in the 
flow and that none of their harmonics or product wavevectors are linearly unstable. These 
secondary waves must be stable even in the presence of finite amplitude primary waves. 

Both expansions are only valid for small amplitude disturbances. Exactly how small is 
not known a priori. In order for the expansions to be valid for large times, it is necessary 
for the linear growth rates to be “small” . Small linear growth rates exist near the neutral 
stability curves. In this sense, the restriction of being near the neutral stability curves with 
the DHZ approach is similar to the requirement in the present approach that the amplitudes 
be small. The present approach might be somewhat more accurate if one were not interested 
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in long-time valid solutions; e.g. if one wanted only to know the first effect of nonlinearity 
on the strength of an instability. In such a case, the present amplitude expansion could be 
valid in cases with large linear growth rates as long as the amplitudes were still sufficiently 
small. 

The theory developed by DHZ describes the weakly nonlinear behavior of the flow for a 
range of Reynolds number and curvature ratio perturbations away from the neutrally stable 
conditions. In particular, the perturbations could be taken into the linearly stable regime; 
there the Herbert approach may fail due to a resonance. The theory of DHZ does not have 
this problem since all of their constants are computed right on the neutral stability curves. 
On the other hand, the present theory may be applied to any combination of wavenumbers 
and Reynolds numbers, while the DHZ theory in its current form works only when the 
azimuthal and axial wavenumbers correspond to neutrally stable disturbances at the same 
Reynolds number. 

The derivation and solution of the equations of the present theory are somewhat simpler 
than those of DHZ since the present theory does not require the use of adjoint eigenfunctions 
and the orthogonality condition. Additionally, because the amplitude function is uniquely 
defined in terms of the fundamental, the change in shape of the fundamental disturbance can 
be computed. This allows the present theory to be extended to higher orders unambiguously. 
Herbert [10] presents seven Landau constants for plane Poiseuille flow using this approach. 
With modified versions of the same procedure used here, Singer, Meyer, and Kleiser [14] and 
Ng, Singer, Henningson, and Alfredsson [15] compute multiple Landau coefficients for the 
development of vortices in a three-dimensional boundary layer and a rotating channel flow, 
respectively. Unfortunately, the derivation of the equations for the higher Landau coefficients 
for the present interacting wave case is too complicated to be handled within the framework 
of our current Macsyma program and is not done here. 


III. Results 


Most of our results are at the critical curvature parameter A c since this is where we expect 
the effects of the interaction to be greatest. Near the critical curvature parameter, A — 
— 2.179 X 10 5 , we found that small changes in the curvature parameter do not change 
the qualitative properties of the results. To orient the reader, Figures 2 and 3 are neutral 
stability curves for the TS and Dean waves respectively at the critical curvature. Regions 
where o 0 , o and 6 0) o are positive represent the linearly unstable areas. All of our results will 
be for points that are in the unstable regime of both the TS and Dean disturbances. 

Figures 4 and 5 illustrate the zero contours of the Landau coefficients ai,o and £>o,i i n 
regions where the TS and Dean disturbances are linearly unstable, respectively. These 
coefficients represent the self-interaction effects of the TS and Dean disturbances. Positive 
values indicate that increased amplitude exasperates the instability. Negative values are 
necessary for an individual disturbance to reach a saturation state in the linearly unstable 
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region. The TS waves have negative self-interaction coefficients near the lower branch of 
the TS neutral stability curve, but positive values elsewhere. At low Reynolds numbers, the 
Dean disturbances have negative self- interaction coefficients for all wavenumbers /3; at higher 

Reynolds numbers the self-interaction coefficient is positive in the middle of the unstable 
region. 

Figures 6 and 7 indicate the effect that interaction has on the disturbances. In both of 
these figures, at a given Reynolds number, the right and left sides of the rectangle bound 
the streamwise wavenumber range in which individual 2D TS waves are linearly unstable. 
Similarly, for each Reynolds number, the top and bottom sides of the rectangle bound the 
spanwise wavenumber range in which Dean disturbances are linearly unstable . In Figure 6, 
positive values of a 0 ,i indicate that the Dean disturbance has a destabilizing effect on the 
TS wave. This occurs primarily in the upper right hand corner of the unstable rectangle, 
though at high Reynolds numbers, there is also a small region of positive a 0>1 in the lower left 
corner. Where a 0( i is positive in the upper right corner, the TS self-interaction coefficient 
a li0 is also positive, so the destabilizing effect of a 0fl here exists when the TS wave by 
itself is already quite unstable. The Dean disturbance generally has a stabilizing effect on 
the TS waves at low Reynolds numbers, but as the Reynolds number increases, the Dean 
vortex has a destabilizing influence over a greater portion of the domain. The corresponding 
results for & 1>0 , the effect of the TS wave on the Dean mode, are illustrated in Figure 7. 
Low wavenumber Dean disturbances tend to be destabilized by high wavenumber TS waves; 
however the fact that the harmonic of the Dean is also unstable in these cases invalidates 
the application of the theory here. In no case does the theory indicate that TS and Dean 
disturbances simultaneously destabilize each other. Since the coefficient, b lfi is positive over 
progressively smaller portions of the domain with increasing Reynolds number, it is unlikely 
that simultaneous destabilization occurs, even at very high Reynolds numbers. In most of the 
cases studied, the weakly nonlinear interaction tends to stabilize both modes. In a physical 
experiment, multiple three-dimensional modes can interact, perhaps forming a secondary 
instability which will eventually grow and lead to transition. A multi-wave interaction and 
subsequent instability is beyond the scope of the present theory. 

For several representative sets of parameters, we examine our results in the context of 
nonlinear autonomous systems. Recall that A measures the amplitude of the TS wave and 
B corresponds to the Dean disturbance. Equations 35 and 36 admit four possible steady 
state solutions: the trivial solution, finite A with 5 = 0, finite B with A = 0, and a 
combined state with finite values of both A and B. Each of these states is a critical point 
and may be classified according to the local behavior of the solutions in the vicinity of the 
point. In the phase plane (A-B), stable nodes are identified as those for which solution 
trajectories near the point converge to the node, unstable nodes are critical points from 
which solution trajectories diverge, and saddle points are those points for which a finite 
number of trajectories converge to the point while all others diverge. 

Several trajectories are illustrated in figure 8 for the case where Re = 6291.67, A = 
2.189 x 10 , & — 74257, and j3 = 4.51. These parameters were used extensively by Singer 
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and Zang [9] in their direct numerical simulations. In the figure, the small squares indicate 
the various initial conditions. An explicit fourth order Runge-Kutta differencing scheme was 
used to trace the subsequent trajectories. The filled circles show the two stable equilibrium 
points at A = 3.8 X 1(T 3 , B = 0, and A = 0, B = 7.0 X 10 -5 respectively. Trajectories 
initiated near either of these points tend to converge to the respective nodes. The open 
circles represent unstable nodes and saddle points. Two initial conditions started quite close 
to the saddle point at A = 2.4 x 10 -4 , B = 3.5 x 10~ 5 go to different stable equilibrium 
states. The flow is especially sensitive to small changes near the saddle point. The unstable 
node at the origin is the endpoint of a semi-infinite curve which goes through the saddle 
point and separates the region of attraction of the two stable nodes. 

Figure 9 illustrates a typical trajectory for the case with Re = 4000, A = 2.179 x 10 , 

a = 96300, and (3 = 1.55. In this case there are three critical points. The origin and A = 0, 
B = 4.69 x 10 -5 are both unstable nodes. Trajectories that start in the vicinity of these 
points move away from them. The point A = 3.8 x 10~ 4 , B = 6.1 x 10 -5 is a stable spiral 
node. Trajectories slowly spiral towards this point. There is no conventional equilibrium 
point with B = 0 and A / 0. Trajectories initiated along this line will tend towards A — » oo. 

In the next case, the flow parameters are Re = 5000, A = 2.179 x 10 -5 , a = 100000, 
and /3 = 2.0. In Figure 10 pairs of trajectories with initial values of A differing by 0.001 are 
plotted for various initial values of B. In many ways, the qualitative behaviour seen here is 
quite similar to that illustrated in Figure 8. The origin is an unstable node, there is a saddle 
point with nonzero A and B, and there is a stable node with nonzero B but A — 0. A semi- 
infinite curve extends from the origin, through the saddle point and separates the domain 
into a region which is attracted to the stable node with A = 0, B = 7.7 x 10 -4 and a region 
in which A — > oo. This shows how strong the nonlinear interaction can be. Without a Dean 
disturbance, the TS wave has unbounded growth; however, the inclusion of an additional 
disturbance at sufficiently large amplitude can completely stabilize the otherwise growing 
TS wave. 

We conclude this section by comparing our Landau coefficients with those of DHZ for 
two specific cases. Tables 2 and 3 show the Landau constants from the current theory 
and from the theory of DHZ after conversion to the current nondimensionalizations. The 
quantities a 0 , 0 and b 0>0 from DHZ were taken directly from their linear eigenvalue solver, 
rather than by perturbing away from the neutral stability curve. Except for the coefficient 
Ji o, which represents the effect of the TS wave on the Dean disturbance, the coefficients 
show satisfactory agreement. Singer and Zang [9] suggested that a yet unknown numerical 
problem in DHZ led to erroneous values of fe 1>0 . The large discrepancy here with respect to 
6 1>0 adds support to that suggestion. 
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a 0,0 

a l,0 

a 0,l 

Vo 

&1,0 

V* 

Current 

1.21 x 10~ 4 

-8.25 

-95000 

3.97 x 10" 5 

-482 

-7640 

DHZ 

1.20 x 10" 4 

-8.85 

-101000 

3.97 x 10 -5 

1.23 x 10 5 

-8470 


Table 2: Comparison of Landau Coefficients for Re = 6291.67, A = 2.189 x 10~ 5 , a = 74257, 
and /3 = 4.508 



a 0,0 

°1,0 

a 0,l 

Vo 

Vo 

Vi 

Current 

1.18 x 10~ 5 

3.97 

-46600 

7.66 x 10" 6 

-233 

-3340 

DHZ 

1.33 x 10" 5 

3.97 

-46900 

7.66 x 10“ 6 

1.44 

-3390 


Table 3: Comparison of Landau Coefficients for Re = 4175, A = 2.179 x 10 -5 , a = 86800, 
and /? - 2.76 


IV. Comparison with Direct Numerical Simulation 


We use direct numerical simulations (DNS) of the full Navier-Stokes equations to determine 
the usefulness of the weakly nonlinear theory. The numerical simulation code [9] uses a curved 
channel variant of the method described by Zang and Hussaini [16] with the nonlinear terms 
in skew-symmetric form [17]. Table 4 reports the complex growth rates (a 0| o — iffo.o) obtained 
with a spectral linear stability code, and the results of DNS with two different time steps. 
Table 5 provides details of the parameters for the different cases. All codes used 65 points 
across the channel. Data from the DNS were obtained after 100 time steps. Using a time 
step At = 0.0001, the complex growth rates differ from those predicted by linear theory by 
less than 1 part in 10 4 . In the simulations reported below, we are primarily concerned with 
long-time trends and final steady-state solutions, hence, the additional time advancement 
errors associated with the larger time step, A t = 0.01, are not expected to be important. 
The larger time step gave a maximum CFL number of 0.01, still well within the numerically 
stable regime. In one case described below, the time step was changed from At = 0.01 to 
At = 0.0001 after the steady state solution was reached. There was no evidence that the 
solution was changing, even after 20000 time steps. Additional details of code validation 
studies in both the linear and nonlinear regimes are given by Singer and Zang [9], 


Case 

Linear Stability 

DNS At = 0.01 

DNS At = 0.0001 

i 

(3.1189,-4236.9) 

(3.0736,-4241.8) 

(3.1185,-4236.9) 

2 

6.6331 

6.6336 

6.6335 

3 

(-138.2126,-16214.84) 

(-138.2154,-16214.84) 

(-138.2126,-16214.84) 

4 

(86.4287,-2709.11) 

(86.4204,-2709.44) 

(86.4291,-2709.14) 

5 

21.5591 

21.5599 

21.5599 

6 

(-90.61188,-12993.29) 

(-90.61307,-12993.3) 

(-90.61188,-12993.34) 


Table 4: Comparison of Complex Growth Rates (all values multiplied by 10 4 ) 
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Case 

Re 

a 

p 

i 

5000 

100000 


2 

5000 



3 

5000 

100000 


4 

10000 

80000 


5 

10000 



6 

10000 




Table 5: Specification of parameters for cases in Table 4 

Sufficient spatial resolution in the simulation can be ensured using the guideline suggested 
by Krist and Zang [18]. They suggest that “grid refinement is needed in any direction when 
the tail of the energy spectrum reaches 10 -8 of the low-frequency value.” This guarantees 
that truncation errors in the velocity will be less than 0.01%. Such detailed resolution is 
not necessary for the purposes of these simulations. We found that the tail of the energy 
spectrum could be as much as 10 -4 of the low frequency value and still provide similar results 
to those obtained when the Krist and Zang guideline was strictly followed. Our relative 
insensitivity to resolution is attributed to the rather minor role that small scales play at this 
stage of transition. In what follows, using M Chebyshev modes implies that there are M + 1 
Chebyshev collocation points across the channel; using N Fourier modes implies that there 
are N collocation points in the given direction so that the Fourier wavenumbers, n, have a 
range of — N/2 + 1 < nL/( 2 -k) < N/2 — 1, where L is the length of the domain. In most 
of the simulations with Re = 6291.67, 48 Chebyshev modes are used to resolve the flow in 
the wall-normal direction, though in the simulations with Re = 5000, 64 modes are used. 
In the streamwise and spanwise directions, 2-12 Fourier modes are used, depending on the 
purpose of the given simulation. The code is structured such that at least 2 modes must 
be used in each direction. In cases where no TS mode is included, 2 modes are used in the 
streamwise direction. Similarly, where no Dean mode is included, 2 modes are used in the 
spanwise direction. 

The initial disturbances for the direct numerical simulations contain only contributions 
from the linear solutions. The distortion of the mean flow and higher harmonic contributions 
are not included initially; they develop as the simulation progresses. 

Direct numerical simulation is used to study two of the three flows described in the 
previous section. The case in which there is a stable spiral mode is not simulated because the 
time required for a single orbit is so large that the computational cost would be prohibitive. 

We first look at the case with Re = 6291.67, A = 2.189 x 10 -5 , a = 74257, and (3 = 4.51. 
Four simulations were performed with different initial conditions; the first contained only 
a Dean vortex, the second contained only a TS wave, the last two contained both types of 
disturbances. 

When only a Dean disturbance is included in the initial conditions, the weakly nonlinear 
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theory suggests that an equilibrium state with B — 7.0 X 10 ® develops. This corresponds 
to a maximum streamwise velocity perturbation of 5.84% of U*. We initiated the simulation 
with the Dean eigenfunction obtained from linear stability theory with an initial strength, 
B = 6.95 x 10~ 5 . Twelve Fourier modes were used in the spanwise direction, though half 
the modes had energies which were more than 8 orders of magnitude less than the primary 
mode. These higher order modes had a negligible effect on the simulation. At the start of 
the simulation, the instantaneous growth rate of the Dean disturbance was approximately 
equal to its linear value. At the time, t = 34845, the magnitude of the instantaneous growth 
rate had decreased by more than a factor of 100 and the simulation was stopped with 
B = 6.7 X 10 -5 . In this case, the difference between the equilibrium state predicted by the 
weakly nonlinear theory and that obtained from the DNS is less that 5% and is considered 
quite good. 

With only a TS wave, we expect to obtain a TS equilibrium state with A = 3.8 x 10 -3 , 
resulting in a maximum streamwise disturbance of 1.94% of U*. The simulation was initiated 
with A = 3.5 x 10~ 3 and was stopped at a time t = 11367. At this time, A = 4.0 x 10~ 3 . 
Figure 11 illustrates the time evolution of the amplitude of the wave. The dashed line 
shows the equilibrium amplitude predicted by our weakly nonlinear theory. Note that as the 
amplitude approaches its equilibrium value, the growth rate decreases. Over 480 periods of 
the TS wave were simulated for the data presented in Figure 11. The difference between 
the equilibrium amplitude predicted by the weakly nonlinear theory and that obtain by the 
DNS is about 10%, which we considered acceptable agreement. The 12 Fourier modes used 
in the streamwise direction were sufficient for an 8 order of magnitude drop in the energy. 

We performed two simulations with both Dean and TS perturbations included in the 
initial conditions. In these, we used 6 modes in both the streamwise and spanwise directions 
and 48 Chebyshev modes in the wall- normal direction. In these simulations, the spanwise 
resolution was only sufficient to allow an energy drop of 10 s , while the streamwise resolution 
allowed as little as a 10 4 drop in the energy. This would lead us to suspect errors in the 
velocity on the order of 1%. To determine whether this could lead to erroneous physical 
results, we re-simulated our case which had the largest velocity perturbations using 12 modes 
in the streamwise and spanwise directions and 64 modes in the wall-normal direction. The 
energy spectrum in the highly resolved case exhibited an energy drop of 10 9 in the streamwise 
direction and 10 14 in the spanwise direction. The evolution of the disturbance amplitudes 
was the same as that observed in the less resolved case. In both cases we initiated the Dean 
disturbances with B = 6.95 X 10 5 . This value is slightly less than its equilibrium amplitude 
which we obtained above. 

In the first case, the initial strength of the TS disturbance was A = 2.94 x 10~ 4 , approx- 
imately 1/13 its predicted equilibrium value. The long time behavior of this case required 
12 Fourier modes in the spanwise direction. The amplitude evolution of the disturbances is 
illustrated in Figure 12a. Here we plot the amplitude of both the Dean and TS disturbances, 
normalized by their respective amplitudes at time t = 0. The Dean disturbance develops to- 
wards its equilibrium state while the TS disturbance decays. In Figure 12b the predictions of 
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the weakly nonlinear theory are plotted. The agreement between the theoretical predictions 
and the DNS results is very good. 

In the second case, the initial strength of the TS disturbance was A = 2.94 x 10" 3 , ap- 
proximately 10/13 its predicted equilibrium value. Figure 13a shows the amplitude histories 
for this case. Unlike the previous situation, here it is the TS wave which grows towards its 
equilibrium value and the Dean disturbance which experiences rapid decay. The predictions 
of our weakly nonlinear theory are plotted in Figure 13b. Again the agreement between 
theory and simulation is quite good. 

The next flow which we investigated with numerical simulations corresponds to the last 
case studied in the previous section. Here Re = 5000, A = A c = 2.179 x 10 5 , a = 100000 
and (3 = 2.0. Simulations were performed with a TS wave only, a Dean vortex only, and a 
combination of both disturbances. 

The weakly nonlinear theory predicts that a TS wave alone in this flow becomes more 
unstable as its amplitude increases since both ao,o an d ai,o are positive. A numerical simu- 
lation initiated with A = 0.002 and B = 0 confirmed that the amplitude A increased with 
increasing growth rate, at least until A = 0.0048 when the simulation was terminated. Six 
Fourier modes were used in the streamwise direction. 

A Dean vortex alone in this flow is predicted to reach an equilibrium state with B = 7.7 x 
10" 4 . This value of B corresponds to a maximum streamwise velocity perturbation which is 
62.4% of the undisturbed laminar bulk velocity. An equilibrium state with B = 4.3 x 10 -4 
was obtained in the direct numerical simulations. Twelve Fourier modes were used in the 
spanwise direction. The difference between the DNS result and that of the weakly nonlinear 
theory is not unexpected since the disturbance amplitudes are so large. 

One final simulation was performed with the initial values of A and B chosen as A = 0.002 
and B = 0.00065. In this case, the weakly nonlinear theory predicts that the TS wave will 
decay, leaving only a Dean disturbance in the flow. The time evolutions of the TS and Dean 
disturbances are plotted in Figure 14a through t = 500. Six Fourier modes were used in 
the streamwise direction, twelve in the spanwise. At this time the spectral energy decay 
in the spanwise modes is 10 5 . Though the simulation was terminated before a steady-state 
solution was reached, the prediction of the weakly nonlinear theory (plotted in Figure 14b is 
qualitatively supported. After the initial transient, the TS wave decays rapidly. Considering 
the large amplitude disturbances in the flow, the theory and the simulation agree well. 


V, Conclusions 


We have presented a weakly nonlinear theory to describe interactions between TS waves and 
Dean vortices in curved channel flow. The approach used to calculate the Landau coefficients 
is an extension of that developed by Herbert [10, 11]. This approach provides not only 
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the Landau coefficients, but also the shape change of the fundamental disturbances with 
amplitude. A standard Chebyshev collocation procedure is used to obtain numerical results. 
We have restricted our study to very weakly curved channels, where the critical Reynolds 
numbers for both TS and Dean disturbances are approximately equal. For low Reynolds 
numbers the interaction effects are generally stabilizing to both types of disturbances. In 
the case of high Reynolds numbers the theory predicts that the interaction can destabilize 
tfi e TS wave. Destabilization of the Dean vortex due to the TS wave occurs only for very 
low wavenumber Dean disturbances. 

Three different scenarios of TS-Dean interactions have been studied in detail. Phase 
plane trajectories indicate the existence of stable and unstable critical points. Two cases 
exhibit saddle point behaviour while the third has a stable spiral node. 

Direct numerical simulations verify the predictions of the theory. Wfien the amplitudes of 
the disturbances are not too large, the direct numerical simulations and the weakly nonlinear 
theory agree quite well. When the amplitudes of the disturbances are large (on the order of 
60% of the centerline velocity), the simulations and the theory still qualitatively agree. 

Comparison of our results with those of DHZ, who also developed a similar weakly non- 
linear theory strongly suggest that DHZ have a numerical error in the computation of one 
of their Landau coefficients. 
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Appendix A 


APPENDICES 


I. The Left Hand Side Operator L n> / 
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II. The Right Hand Side Vectors R s 


R, = 


/*i\ 

R] 

Rl 

\rV 


by: 


Equation system 5 corresponds to the TS harmonic. Its right hand side vector is given 


R 5 


R\ 


R\ 


Rt 


— ia ug 1 ,0,0 ’ 0 u r 1 ,0 ’ 0,c y r 
+(u < , 1 - 0 ’ 0 ' 0 ) 2 /r 

-ia(u e lfi ’ 0 '°) 2 /r 
— Ug 1 ’ 0 , 0 , 0 n r 1 ’ 0 ’ 0,0 / r 


(Ug 1,0 ' 0 ' 0 )u T 1 ’ 0,0 '° 

= -iaug^^u^/r 


dr 


+UT 1 ' 0 ' 0 ' 0 ^(u, 1 ’ 0 ’ 0 ’ 0 ) 

dr 


= 0 


by: 


Equation system 6 corresponds to the Dean harmonic. Its right hand side vector is given 


R\ 


= —ifiu °>°> 1 ' 0 - °' 0 ’ 1 - 0 
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+(u e °’ 0 ' 1 '°) 2 /r 

dr 
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dr 


R 6 


Rt 


= -i/^ 0 ' 0 ' 1 - 0 ) 2 
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0,0, 1,0 


dr 


a 0,0,1,0\ 


u 


= 0 


Equation system 7 corresponds to the product of the TS and Dean fundamentals. Its 
right hand side vector is given by: 


R\ ~ -iau e 0 i ° tl, 0 u T 1,0,0,0 /r 
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-i/K°'°' 1 ,( V i0 ’ 0,0 

+2u e 0 ’ 0 ’ 1 ’W fi ' 0 '°/r 
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R$ = 0 


Equation system 8 corresponds to the product of the TS fundamental and the complex 
conjugate of the Dean fundamental. Its right hand side vector is given by: 
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+z/3u r 0,0,-1,0 u x 1, 0,0,0 
+ 2 t 2 fl 0 , 0 ,- 1 , 0 ttfl 1,0 ’ 0,0 /r 

- ( -7-(u r 0,0,_1 ,0 )u r 1 ,0, °’° + ^-(u r 1,0,0,0 )zz r 0,0, - 1, °) 
dr dr 

« - 0 , 0 , - 1,0 A 1 , 0 , 0,0 # 

— zaiifl ’ ’ ' ue /r 

0 , 0 ,- 1,0 ; 1 , 0 , 0,0 

Ti_ 


-(tx fl 0,0 '- 1,0 u r 1, 0,0,0 + h fl 1,0,0,0 u r 0,0, - 1,0 )/r 
v dr dr 


• - 0,0, -1,0 * 1,0,0, o / 

—iau e ’ ’ ’ u z /r 


7 oa-i.o^ l, 0 , 0,0 


-f z/3u 

_(A(^o.o,-i.o )tIr i, 0,0,0 + A( U V ,0, 0 , >V ), 0 , -1’0) 

v dr dr 


= 0 


Equation system 9 corresponds to the self correction of the TS fundamental. Its right 
hand side vector is given by: 

Rl = — za(2uV 1,0 ’ 0,t V l0 ' 0, ° + U0 0 " 1, 0,0 u, 1, 0,0,0 - txfl 2,0,0,0 t2 r _1,0,0,0 )/r 


+ (2ur 1 ’ 0 ' 0 ' 0 ^ 2,0,0,0 + 2 u (? °’ 1 ' 0,0 ^ 1,0,0,0 )/r 


_( A( u ;2.W )l2r -i.o,o,o + “ (l2p -i.o.o.o )tip2 
v dr dr 


-i,o,o l o\, ; 2 , 0 , 0,0 
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Rl 


+ ^.( l 2 p O,i 1 o,ow 1,0,0, o + A( t 2 r 1 - 0 - 0 -0) t i r 0 ' 1 -o-0) 

dr dr 

-K 0 - ^i.o)^ 1 ’ 0 ' 0 ' 0 


JK = 


-ia(u 9 - 1,0 ’ 0, % 2,0,0 ' 0 + u* 

°u T 2 


0.1.0, 0-1 ,0,0,0 


)/r 


U0 

-( nV 1 ' 0 ' 0 ' 0 ^ 2 ' 0 ' 0 ' 0 + liV ' 0 ' 0 ' 0 ^" 1 ’ 0 ' 0 ’ 0 + us 0 ' 1,0,0 ^ 1 ' 0,0,0 + Ufi 1 , 0 ' 0 , 0 ti r 0,1 ’ 0,0 )/r 


+ it*-— -tx, 

_(i.(^-i.o.o.o )t2r 2, 0,0,0 + 
ar dr 

+ 4-(ue 0 ’ 1 ’ 0 ' 0 )u r 1 ’°' 0>0 + ^(uV' 0 ’ 0 ’ 0 )^ 0 ’ 1 ’ 0 ’ 0 ) 
ar ar 

-(ai.o - ^i,o)uV' 0 ’ 0 ' 0 


= —ia(2ug~ 1,0,0 '°u z 2 ' 0 ' 0 ' 0 + u e 


0 , 1 , 0,0 - 1 , 0 , 0,0 _ « _ 2 , 0 , 0,0 * - 1 , 0 , 0,0 


u 


— u e 


U r 


3 )/l 


_(.i(u z - 1 . 0 . 0.0 V 2 . 0 . 0.0 + iu/MOw-iMO 

ar ar 


d 


d 


+ ^(^0,1,0,0)^1,0,0,0 + ( 

dr dr 

-Ko - ^.o)^ 1 - 0 * 0 - 0 


- 1,0, 0,0 N . 0,1 ,0,0 


ix 


K 0 ’ 1 ’ 0 ’ 0 ) 


RX = 0 


Equation system 10 corresponds to the self correction of the Dean fundamental. Its right 
hand side vector is given by: 

R\ 0 = —i/3(2u z 0,0, ~ 1 ' 0 u r 0,0,2 '° + u z 0 '°- 0 - 1 u r 0 - 0 - 1 - 0 - u z °' 0 - 2 - 0 u r 0 ' 0 --i-0) 


Rio 


R 3 


10 


-\-{2ug 


0 , 0 , - 1 , 0 ^ 0 , 0 , 2,0 


+ 2tt e 0,0 ’°’ 1 ■u fl 0,0,1, °) / 


ar ar 

+ A(u r 0.o.o.i)u r °.o.i.o + A(u r 0 -°. 1 - 0 )u r 0'0.0,i) 

dr dr 

-(6o,i - i/>o,iK 0,0,1, ° 

— i/3(2i2fl 0,0 ' 2,0 u z 0,0 ’ _:l ' 0 + ix, 0 - 0 - 1 - 0 ^ 0 - 0 - 0 ’ 1 - ^o.o.-i.o^o.o^.O) 

- 1 > V 0,0 ' 2 ’° + tx e 0 - 0 - 2 - V 0,0, “ 1,0 

+u7- 0A1 u r °' (U - 0 + u e 0 ' 0 ' 1,0 u r 0,0,0 'i)/r 

_(£(^ 0 . 0 .- 1 . 0 ) tir 0 , 0 . 2,0 + A (^ 0 , 0 , 2 . 0 ^ 0 , 0 , - 1,0 

ar ar 

+ i.(^o.o,o,i )t2r 0 , 0 , 1,0 + d o,o,i,o )tIr O,o,o,i ) 
dr ar 

-(Vi - ;vi)«« 0,0,1, ° 


-ip(u z 0fi '- l ' 0 u x 0 ' 0 ' 2 ' 0 


0,0, 0,1 y 0,0, 1,0 


+ u z -’-’-’- U 

d 

dr 
d 


) 


7 0, 0,-1, 0)^0, 0,2,0 _|_ _^_^0, 0,2,0)^ 0,0, -1,0 


+ ^(O* 0 -i)u r 0 * 0,1,0 + ^o.o.i.o^o, 0,0,1) 
ar ar 

— (^0,1 — iho>)u z °W 
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0 


Rt - 


*-10 


Equation system 11 corresponds to the correction of the TS fundamental by the Dean. 
Its right hand side vector is given by: 

= -iaiue 0 ' 0 ’- 1 ' 0 ^ 1 ' 0 ' 1 ' 0 + ue 0 ’ 0 ' 0 ' 1 ^ 1 ' 0 ' 0 ’ 0 + u, 0 , 0 , 1 , V 1,0, - 1, °)/’- 

+i/3(u r 0,0,_1,0 u x 1 ’ 0 ’ 1 ’ 0 - tir 0.o.i.o tij 1 .o.- 1 .° 

+tt P 1, ° ,-1,0 '“* 0 ’°’ 1, ° ~ iI r ll °* 1,0 O 0,-1, °) 

+2(u^°'- 1, % 1, 0,1,0 + u* 0 ’ 0 ’ 0 ’ 1 ^ 0 ’ 0 - 0 + ug O,o, 1 , °U 0 L,o, ~ 1,O )/‘r 


-^U^'-^Ur 1 ' 0 ' 1 ’ 0 + ^(Ur 0 ' 0 ' 0 ' 1 )^ 1 ' 0 ' 0 ' 0 

x dv &T 

dr dr 

+ A( a 2 r 1 ' 0 ' 0 '°)u r 0 ' 0 ' 0 ’ 1 + ^(u, 1 ' 0 - 1 - 0 )^ 0 ' 0 -- 1 ’ 0 ) 

dr dr 

— (®o,i - ^o.i )^ 1 ’ 0,0 ’ 0 

R\ x = -iaCus 0 ’ 0 ’- 1 ’ 0 ^ 1 ' 0 ' 1 ’ 0 + ue 0 ' 0 ' 0 ' 1 ^ 1 ' 0 ’ 0 ’ 0 + u 6 ^W’ 0 '-''°)/r 


+ Ug^'^Ug^’^ + Ug Ug' 

+ip(ug 0,0, ~ 1 ’ 0 u z 1,0 ' l ’° - ue 0fiX0 nV* 0 -- 1 ’ 0 


+u fl 1 ’ 0 '- 1 ' 0 ^ 0 - 0 ' 1 ' 0 - ue 1 ' 0,1,0 ^ 0,0 ’ -1,0 ) 


■ 0 > 0 '- li0 t£ r 1 * 0>1<0 + ue 0,0,0,1 ^ 1, 0,0,0 
r _0,0,l,0, ? 1,0, -1,0 + ^l.O.-l.O^O.O.l.O 


~{ug 

+ Ug 

+ Ug'' 0 '°'V' 0 ' 0 ' 1 + Ug 


U T 


• 1 ’ 0,1 ' 0 i2 r 0,0,_1 ’° 


*)/r 


— ( — (u 9 0,0 ’- 1,0 )u r 1,0,1,0 + A( u V 0 - 0, l)^r 1,0l0 ’0 

v dr dr 

+ ±U 2 j? 0,0 . 1 .°)u r 1 ' 0 '- 1 - 0 + i-( U V ,0, " 1,0 )«r 0, 0,1,0 

dr dr 

+ A( t 2 fl 1 .°.o,o )tir 0.o.o,i + ^.( u V 0ll ’ 0 )^ 0A " 1,0 ) 
dr dr 

-(a 0 ,i - igorf * 1 * ),0, ° 

R 3 n = —ia(u0° ,o, ” 1,o u* 1, 0,1,0 + ufl 0,0,0,1 ^* 1, 0,0,0 + u e 0, 0,1,0 «x 1,0, ' 1,0 )/r 

-(—(u°' 0 ~ 1 ‘ 0 )Ur 1,0,1 '° + ^-(u* 0 * 0 - 0 - 1 )^ 1, 0,0,0 
v dr dr 

+ 1(„- 0.0,,.0) Ur 1.0,- 1 .0 + i 4 «-W)4WI 

dr dr 

+ A(uV.°.0.0w 0.0.0,! + ^.(^1.0,1. 0)^0, 0,- 1( 0) 

dr 


/ * \ -7 1 , 0 , 0,0 

— ( a o,i - * 2 o,iK 


itfi = 0 


Equation system 12 corresponds to the correction of the Dean fundamental by the TS. 
Its right hand side vector is given by: 

R' u = -ia(u e -' A0 -°Ur l - W + ir 1 "* 1 *" - u, 1 ' 0 ' 0 ’ V 1 "' 0 + u, IM '°u r -'' W )/r 
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R 12 


*?2 




-i(3{u r W-%1.0,0,0 + ^0,0,1, 0^0,1, 0,0 + u r 1.0.1, 0 U —1, 0,0,0) 

+2(t2r 1,0,0 *% 1>0>1>0 + uV 1 ’ 0 ' 1,0 ^ 1 ’ 0 ' 0 ' 0 + uV’ 0 ' 1 ’ 0 u e 0 ’ 1,0,0 )/r 

-(^r(^r' 1,0,0,0 )^ 1,0,1,0 + ^-(t/ r - 1 > 0 ’ 1 - 0 )u r 1 - 0 ' 0 > 0 

or 

+^K 0,0,1,0 K 0,li0 '° + -^-(ti P 0,1,0,0 )t2 r 0,0,1, ° 

fl/ 

+ ^^“ rl ' 0,0 ’^“ r ~ 1,0,1 ' 0 + jr(Wr 1,0 ’ 1,0 )Wr" 1,0 ' 0, °) 

-(^i,o - ^i,o)u r 0,0,1 - 0 

= — i / 3 ( Ufi ~ 1 , 0 , 1 , 0 ii 2 1 ,0,0, ° + u /' 0 ' 1 - 0 ^ 0 - 1 ' 0 - 0 + - 1 . 0 . 0 . 0 X 

-(« r 1 ' 0 ’ 0 - V 1 - 0 ' 1 - 0 + ur 1 ' 0 ’ 1 - 0 ^ 1 ’ 0 ’ 0 - 0 

+U0 O '°' 1 '°u r 0, 1,0,0 + uV ),1,0,0 u r 0, 0,1,0 

W'° A V ~ 1>0>1 ' 0 + -i.o.o.O) /r 

-(^(u*~ 1,0 ’ 0,0 )Ur 1 ‘ 0 ' 1 '° + ^ r (^- 1 '°' 1,0 )u r 1 - 0 - 0 ' 0 

+ ^r(^ 0 ' 0 ' 1,0 )«r 0,1 ' 0 ' 0 + i-( t 2 9 °.l,0,0^ 0, 0 , l,0 

+ Jr (^ 1 ’ 0 ’ 0 , 0 )« r " 1 ’ 0 ' 1,0 + ^( uV - 0 - 1 ’ 0 ).;- 1 .^ 0 ) 

-(ii,o - i/ii,o)^e 0,0,1 ’ 0 

= —iaCug -1 ’ 0,0,0 ^ 1, 0,1,0 + U V 1 > 0 ll> V >0l0 -° - ^ 1 , 0 , 0 , V - 1, 0,1,0 _ u fl 1,0,1,0 u - 1,0,0,0 Vr 

-ip{u z 1,0,1 > °u z 1 ' 0 ' 0,0 + U Z 0,0,1,0 U % 0, 1,0,0 + ^1,0, 1,0^-!, 0,0,0) 

_ ^^ rl,0,0,0 ^ Arl,0,1, ° + ^(“z" 1,0,1 ' 0 )^ 1, 0,0,0 

+ ^(^ 0,0,1,0 )^ 0llA ° + ^(n, 0,1,0,0 )^ 0, 0,1,0 

+ ^(u z 1 ’° A0 u r - 1 ’ 0 ’ 1 ’ 0 + A( u V 0 -i-0)u r - 1,0,0, 0) 

~(Ko - *Ai,oK°- 0,1 ’° 

= 0 
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Figure 1: Curved channel flow geometry. 
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Figure 6 



80000 


Figure 6: Zero contour of a 0t i, the effect of the Dean disturbance on the TS disturbance. The 
rectangles represent the range of a and (3 over which individual TS and Dean disturbances 
are linearly unstable at the given Reynolds number. 
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Figure 7 



a 

Figure 7: Zero contour of ^.o, the effect of the TS disturbance on the Dean disturbance. The 
rectangles represent the range of a and j3 over which individual TS and Dean disturbances 
are linearly unstable at the given Reynolds number. 
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Figure 10 
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Figure 10: Trajectories in phase-space with Re = 5000, A = 2.179 x 10 -5 , a = 100000, and 
P = 2.0. Squares represent initial conditions, filled circles represent stable critical points. 
Open circles represent other critical points. 
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